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Abstract 

To investigate the consequences of component confinement such as at a glass transition and the 
well-known energy or enthalpy gap (between the glass and the perfect crystal at absolute zero, see 
text), we follow our previous approach [Phys. Rev. E 81, 051130 (2010)] of using the second law 
applied to an isolated system Sq consisting of the homogeneous system S and the medium S. We 
establish on general grounds the continuity of the Gibbs free energy G{t) of S as a function of 
time at fixed temperature and pressure of the medium. It immediately follows from this and the 
observed continuity of the enthalpy during component confinement that the entropy S of the open 
system S must remain continuous during a component confinement such as at a glass transition. 
We use these continuity properties and the recently developed non-equilibrium thermodynamics 
to formulate thermodynamic principles of additivity, reproducibility, continuity and stability that 
must also apply to non-equilibrium systems in internal equilibrium. We find that the irreversibility 
during a glass transition only justifies the residual entropy 5r to be at least as much as that 
determined by disregarding the irreversibility, a common practice in the field. This justifies a non- 
zero residual entropy S'r in glasses, which is also in accordance with the energy or enthalpy gap 
at absolute zero. We develop a statistical formulation of the entropy of a non-equilibrium system, 
which results in the continuity of entropy during component confinement in accordance with the 
second law and sheds light on the mystery behind the residual entropy, which is consistent with 
the recent conclusion [Symmetry 2, 1201 (2010)] drawn by us. 
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I. INTRODUCTION 



A. Irreversibility 

In the previous two papers , to be referred to as I and II respectively here, we 
initiated a general statistical mechanical investigation of non-equilibrium systems undergoing 
irreversible changes during approach to their equilibrium state. It is a very general approach 
that is not restricted to glasses alone but covers all non-equilibrium situations. In I, we 
deal with homogeneous systems, while in II we deal with inhomogeneous systems. The 
approach is truly statistical mechanical in nature, and is based on applying the second law 
of thermodynamics to an extremely large isolated system, which we denote by Sq; it consists 
of the macroscopic system S of interest in a medium denoted by S containing it; see Fig. [H 
We will use body in this work to refer to any of the three systems. According to the second 
law, the entropy Sq of an isolated system Eq, which is a sum of the entropies S and S of the 
system and the medium, respectively, 

Soit) = Sit) + Sit), (1) 

can never decrease in time 3|-[8|: 

^ > 0. (2) 

dt - ^ ^ 

Any change in the entropy of the isolated system is due to irrevesibility, so that d[SQ = 



HI]. Even for a body, which is not 



dSo > 0, and deSo = in the standard notation 
isolated, d^S = dS — dcS > 0, with d^So ^ 0; indeed, the latter can be of any sign. Two 
important observations about the entropy 5*0 are in order for Eq. ([2]) to have a content.: 

1. The entropy 5*0 exists under all conditions; thus it includes the situation when Eq is 
not in equilibrium or when S undergoes a glass transition. 

2. The entropy 5*0 is a continuous function of all its arguments Zo(t) (see below), not 
shown in Eq. ([2]), and the time t, which is shown. 

What happens inside the isolated system (loss of ergodicity in parts of Eq due to compo- 
nent confinement, chemical reactions, phase changes, turbulence, viscous deformation inside 



it, etc.) cannot affect the direction of the inequa! 
pie of non-equilibrium thermodynamics 



ity, which makes it the most general princi- 



All the above entropies, in addition to 
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FIG. 1: Schematic representation of a system T, and the medium S surrounding it to form an 
isolated system Sq. The medium is described by its fields To,Po, etc. while the system, if in 
internal equilibrium (see text) is characterized by T{t), P{t), etc. 



being an explicit function of time t, are also functions of extensive observables and internal 
variables (see II for more details) that themselves may vary with time. Of the observables, 
the energy, volume and the number of particles are the most important extensive observ- 
ables from an experimental point of view. In the sequel to the above two paper I and II, we 
focus on non-equilibrium thermodynamic principles dictating the behavior of a homogeneous 
system S, the general way its non-equilibrium entropy should behave during relaxation, the 
way S is determined in terms of its microstates, and the wa y S must behave when the system 



becomes confined to a component of the phase space |12|-|l4| during some interval of time 
(for example, at a glass transition), and the concept of the residual entropy, the entropy 
that S will have at absolute zero. As our interest is in investigating the out of equilibrium 
properties of S, we will assume throughout this work and as was done in I and II that 
the medium is in internal equilibrium at all times. As a consequence, S has a well-defined 
constant temperature Tq, pressure Pq and other fields and affinities (see below) at all times. 
In contrast, the system may have no well-defined temperature, pressure, etc. unless it is 
also in internal equilibrium, in which case, we will denote their instantaneous values by 
T{t), P{t), etc.; see I for details. They are usually different from the those of the medium 
unless equilibrium is reached, which occurs as t — oo. It was shown in I that T(t) > Tq in 
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a cooling and T(t) < Tq in a heating experiment, and that 

dQ = T{t)dS{t) = TodeS{t), (3) 

where deS{t) represents the entropy exchange with the medium; the irreversible entropy 
generation or the uncompensated transformation of Clausius and de Bonder IJInl within 
the system is represented by diS(t) = dS(t) — dgSit) > 0. 



B. Determination of Entropy 

The entropy can now be calculated in principle by using Eq. ([3]) even when the process is 
irreversible. As the experimentalist has a control over Tq,Pq, the entropy can be calculated 
by varying Tq, keeping Pq fixed. The heat capacity Cp at constant Pq in terms of the 
enthalpy H{t) and the volume V{t) is given by 

Cp(t) 



see I, so that 



Cp < CpdTn/Tn in cooling 

dS{t) = ^dTo - ^ °/ ° \ (5) 

[ > CpdTo/To in heating 
and can be used to obtain the entropy using the standard integration process 

< 8(0) + f CpdTo/To + ASc in cooling 

^0 ' (6) 







> S{0) + f CpdTo/To + ASc in heating 







where 5(0) is the entropy at absolute zero, and A^c represents the sum of all possible 
discontinuities over the range (0,To) mainly due to latent heats at possible first-order tran- 
sitions; it also includes any possible entropy discontinuity due to component confinement. 



The above inequalities are a simple extension of is a well-known classic result [15|, which 
has been rederived in our approach. Because of a non-negative heat capacity, the entropy is 
a monotonic increasing function of the temperature: 

for all finite temperatures. The presence of an internal order parameter can also be accounted 
for as we will discuss later. We find that the integral 

JcpdTo/To > (8) 





provides an upper bound to the entropy S(t) of the system during a coohng experiment, and 
a lower bound for it during a heating experiment, when the system is out of equihbrium. The 
difference in the two bounds, which we will denote by ASb(Tq, Pq) here, strongly depends on 
the rates of cooling and heating. The two bounds are exactly identical so that ASb{Tq, Pq) = 
0, when the system is in equilibrium. Therefore, ASb(Tq, Pq) can be used as a measure of 
the irreversibility when no internal variables (order parameters) are invoked, as was recently 
done by Johari and Khouri [l6||. However, introducing internal variables allows us to obtain 
the entropy without any error, as we will discuss here, even when the system is out of 
equilibrium and irreversibility is present so that the thermodynamic entropy can be always 
calculated, at least in principle. Irreversible entropy generation appears through the internal 
order parameter in the calculation. This is a well-established approach to non-equilibrium 
thermodynamics jsUl; see also [l, 2|, which we carefully discuss in Sect. IIIII for the sake 
of continuity. 



C. Non-equilibrium Thermodynamic Principles 



We have a through understanding of thermodynamic principles (additivity, reproducibil- 
ity, uniqueness and stability, all defined later) and the statistical entropy for an equilibrium 
system |3|-[7|. The situation is not so obvious when the system is out of equilibrium. We 
wish to examine the conditions necessary for the applicability of these principles to, and ob- 
tain the statistical formulation of the entropy of, a non-equilibrium system with the aim to 
develop a non-equilibrium statistical thermodynamics. While in traditional non-equilibrium 
thermodynamics loillll. the entropy of any body is postulated to exist under the as- 
sumption of local equilibrium, its actual value and its statistical connection with microstate 
probabilities is of no interest as only the difference in the entropy has any meaning. The 
same is also true of the thermodynamics applied to glasses jl^ . [17I2OI. The actual value of 
S is normally set by iraposing additional requirement such as by invoking Nernst postulate 



of the third law 



21 



22| . Our interest here is different. We wish to identify the entropy of 



a body when it is out of equilibrium using statistical concepts. It should be stated here that 
he entropy for a body in equilibrium is well understood in equilibrium statistical mechanics 
3|-|6|. Thus, any new attempt to define the entropy for a non-equilibrium state must reduce 

the the well-known form of equilibrium entropy. 
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D. Component confinement, Entropy Reduction and the Residual Entropy 



1. Component Confinement and the Gap 



Phase space confinement at a phase transition such as a liquid-gas transition is a well- 



23| . The component con- 



known phenomenon in equilibrium statistical mechanics ji, 5, 8, 
finement also occurs when the system undergoes symmetry breaking such as during magnetic 
transitions, crystallizations, etc. In such cases, the confinement occurs under equuilibrium 



conditions, and is well understood. The confinement in phase space can also occur under 
non-equilibrium conditions, when the observational time scale Tnhs becomes shorter than 
the equilibration time of the system, such as for glasses whose behavior and 



properties have been extensively studied; see Fig. |2l which is adapted from Fig. 10.2 on 
p. 442 in IJ]. The behavior of the entropies and their consequences for a given cooling 
rate r have been discussed by us in Ch. 10 in [l^, and is based on the energy gap model 
for glasses discussed there. The form of all these entropies is to ensure that they remain 
concave functions of the energy. The basic premise of the model is the well-known fact {4] 
that there exists a gap (an energy gap) between the energies of the perfect crystal and the 
glass at absolute zero, which is at least as much as i?K; it is equal to -Ek for the ideal glass. 
The gap is due to defects that are present in the glass with respect to the perfect crystal 



that raise the energy or the enthalpy of the former |14j . |24| . The additional energy of a glass 
_Ene — -E'k in a glass, however, is due to the presence of excitations relative to the ideal glass. 
It should be observed that the slope of FG at G is infinite to ensure that the energy of the 
system described by FG is -Ene > -^k at absolute zero. The investigation of the energy 
gap model in [l4| is very relevant here as it is this non-equilibrium component confinement 
during the glass transition that we wish to study in this work. The entropy curve FG in 
Fig. |2] denotes the entropy of some arbitrary time- dependent non-equilibrium state includ- 
ing glasses, and is shown disconnected from the entropy curve BAO describing the entropy 
of the disordered state which includes the liquid and the supercooled liquid. The latter is a 
stationary metastable state below the melting temperature Tm, with the crystal representing 
the stable state below Tm. The entropy of the ordered state, which results in the crystal 
below Tm, is schematically shown by BAO. It is a common practice to call the supercooled 
liquid as an equilibrium state when the interest is mainly to study non-equilibrium states 
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FIG. 2: Schematic form of the communal entropy DACK of the disordered hquid/supercooled 
hquid and the entropy BAO of the corresponding ordered crystalhne state as a function of the 
energy; other extensice observables are held foxed. The common tangent to DACK and BAO gives 
the inverse melting temperature Tm- The entropy of some non-equilibrium disordered state for 
some non-zero cooling rate r is given by FG, which must strictly lie below in accordance with the 
second law. It terminates at G at an energy E'ne, which is strictly higher than the energy (taken 
to be here) of the pure crystal at absolute zero or the energy E'k at K on DACK. The entropy 
at G gives the residual entropy 5r. Even though we have shown FG disconnected from DACK for 
some arbitrary non-equilibrium state, it continuously emerges out of the latter at some point C at 
a temperature below Tm in a glass transition. The portion CG in this case represents the entropy of 
the glassy state. As the cooling rate becomes slower, CG moves towards CK until it finally merges 
with it as r — > 0. The terminal point K is known as the ideal glass transition where the entropy 
of the equilibrium supercooled liquid vanishes. As FG moves towards ACK, i?NE moves towards 
the ideal glass energy i?K and the residual entropy continues to decrease until the latter finally 
vanishes in the ideal glass transition at K. However, for any r > 0, S'r remains non-zero due to the 
presence of excitations in the ideal glass that raises its energy above i?K- 
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denoted by FG; the terminology is obviously incorrect if we need to bring the crystal into the 
discussion. As our interest in this work is to mainly study time- dependent non-equilibrium 
states, we will invariably call the supercooled liquid as the equilibrium liquid, knowing well 
its limitation. In a glass transition from the supercooled liquid, the curve FG emerges out of 
BACK continuously at some point C, which is determined by the experimental setup such 
as the cooling rate r. 

As the system is not allowed to probe all the disjoint components in time, it can be 
investigated by our recently developed non-equilibrium approach in I and II; see for example, 
Eqs. dSM])- One of the important observations of such systems is the existence of a residual 
entropy S'r, such as the entropy at the point G in Fig. [2J It is the entropy a glass or a crystal 
would have at absolute zero, and is related to the possible number of disjoint components 
in the phase space [s], Q, 13]. It is merely a measure of the number of lowest energy 



microstates of energy E^-^ > for a given mode of preparation at absolute zero. The 
macrostate of the glass or the crystal is a collection of all these microstates along with their 
probability of occurrence for a given mode; see ^, 23]. The mode of preparation determines 
the energies of these microstates, which need not be the lowest possible energy (such as E^ 
for the supercooled liquid if we are concerned with a glass) for the system IJ]. The lowest 
possible energy microstates represent the equilibrium macrostate at absolute zero. It should 
be emphasized that there is no reason for the lowest energy macrostate to be non-degenerate 
in principle, although it seems to be true for a majority of the systems. Therefore, we will 
assume in this work that Seq,K = for the supercooled liquid at K. This is borne out of 
several exact calculations as discussed elsewhere The system (glass or crystal) can be 
in any of lowest energy microstates (and not the lowest possible enrgy microstates) for a 
given mode of preparation, and the residual entropy is determined by all these microstates 



2. Residual Entropy 

The existence of a non-zero residual entropy (5^ > 0) is very common in Nature, and 
does not violate Nernst's postulate, as the latter is applicable only to equilibrium states with 
a non-degenerate ground state Sect. 64]. Its existence was first demonstrated by Pauling 



and Tolman 25|; see also Tolman [26']. In addition, the existence of the residual entropy 
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has been demonstrated rigorously for a very general spin model by Chow and Wu 27|. For 



ice, the residual entropy was observed by Giauque and Ashley 



Pauling 



29| provided 



the first numerical estimate for the residual entropy for ice, which was later improved by 



Nagle 



30|. Nagle's numerical estimate has been recently verified by simulation 
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32|. The 



numerical simulation carried out by Bowles and Speedy [33| for glassy dimers also supports 
the existence of a residual entropy. Thus, it appears that the support in favor of the residual 
entropy is quite strong. We wish to emphasize that what is customarily called the third 
law due to Nernst, according to which the entropy must vanish at absolute zero, is merely 



a postulate and not a strict theorem even in equilibrium [4, |21|, |22|. Indeed, many exactly 
solved statistical mechanical models show a non-zero entropy at absolute zero. Based on 
these empirical observations, the residual entropy 5*^ at G in Fig. |5] has been taken as 
non-zero and positive. 

We refer to Fig. [2] to emphasize this point once more. As r decreases, G moves not 
only to the left but also down in the figure. This has the consequence of lowering i?NE and 
S'r, so that entropy continues to decrease with temperature to ensure stability of the glass. 
The additional energy or enthalpy due to excitation compared to the ideal glass raises the 
residual entropy of the glass relative to the ideal glass: 



Sr > SenK = for Ene > Ek- 



(9) 



However, it is possible that the residual entropy Sr at G remains zero for all r as if 



all glasses are identical 
[2] has been challenged 



at absolute zero for any r. Indeed, the situation shown in Fig. 



3J-[38|, because it is argued that the entropy cannot be estimated 



correctly in the glass transition region, the region where component confinement occurs, 
where irreversibility comes into play. This is certainly true, as can be easily seen form Eqs 
([6]) and (1601) . While the idea that the irreversibility raises some concern about the inferred 
values of the entropy is certainly justified, the main point is not whether irreversibility is 
present; see Sect. IVI F[ Rather, the question should be: how much of an error does this 
create in the inferred values of the entropy? Thus, the issue at hand is the magnitude 
of AS'sfTn, -Po)- Various estimates support an error that does not exceed more than 5%. 
Jackie ISj has also calculated the amount of uncertainty in the entropy for glycerol due 
to irreversibility near the glass transition, and found it to be less than 5%. Gutzow and 
Schmelzer 39| have shown by citing some earlier analysis that the irreversible contribution to 
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the entropy in the vicinity of the glass transition is small enough to be practically neglected. 
They also provide a very enlightening historical review of the residual entropy in glasses. 



Nemilov 



40| comes to a similar conclusion by suggesting that the error in S'r is no more than 



5%. The conclusion is that there is no appreciable error in using conventional equilibrium 
thermodynamics by replacing the inequality to an equality in Eqs. ([MS]) in studying glasses, 
a point also made recently by Goldstein |4l|. In particular, one can get a highly reliable 
estimate of the entropy of a glass (GL) or the supercooled liquid (SCL) without using any 
internal order parameter by the following standard integration process using the measured 
heat capacity 

S{To) = SiO) + J^^ + ASc. (10) 



More recently, Johari and Khouri 16(| have come to a similar conclusion by reanalyzing 
experimental data. 



3. Entropy Reduction 



It has been suggested |34l438l| that as a consequence of the component confinnment, the 
system is no longer able to explore in time the entire part of the phase space that would be 
consistent with the observables (such as the energy) of the system. With this interpretation 
in mind, and recognizing the fact that the entropy is measured by the volume of the accessible 
phase space, such as that belonging to a component, it is then argued that the entropy S 
of the system must undergo a sudden or rapid reduction in its value due to component 
confinement 138:] that is expected to occur at the glass transition, but the experiments are 
unable to capture this because of the irreversibility. We will call this entropy reduction 
scenario the unconventional view (UV) of the glass transition in this work. In this scenario, 
either FG in Fig. |2] remains disconnected from BAG to account for a sudden discontinuous 
reduction or touches BAG at some point C but falls off rapidly from C to give rise to a rapid 
but continuous drop. In the latter scenario, the portion of FG to the immediate left of C 
will have to change its curvature so that it will be no longer concave, a point of critical 
importance that will be discussed at several occasions in this work. It is impossible to 
continuously join two disjoint concave functions without creating a loss of concavity. 

If it happens that the entropy S of the medium concurrently undergoes a sudden or rapid 
rise by an amount such that it not only compensates the reduction in S, but add some more 
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so that 5*0 does not decrease, then the rapid reduction in 5* would not violate the second 
law for the isolated system as given in Eq. ([2]). However, in this case also, the irreversible 
entropy change in the system must still satisfy the second law: diS{t) > 0. Therefore, our 
third aim is to investigate if such a scenario can be supported by the second law. 

Our last goal is to use the second law to investigate the idea of the residual entropy 
from a statistical point of view and to prove that the situation in Fig. [2] with a non-zero 
residual entropy is borne out of experiments. This will be in contradiction with the idea of 
an entropy reduction in the unconventional view. The enthalpy H and the volume V are 



known to be continuous during the glass transition Thus, any entropy reduction 

will leave its imprint on the behavior of the Gibbs free energy. Thus, we will also investigate 
the behavior of the Gibbs free energy in this work. 

The other relevant issue, which to the best of our knowledge has not been discussed 
in the literature so far for unexplained reasons, is whether the experimental data support 
a non-zero residual entropy even if the irreversibility may be quite appreciable, as would 
happen in a rapid quench. From above, it is clear that the experimental determination 
of the residual entropy S'expt(O) invariably yields a non-zero value. We derive later a very 
simple consequence of the glassy irreversibility in the form of an inequality; see Eq. fl60|) . The 
inequality clearly supports a non-zero residual entropy, even if the irreversible contribution 
is appreciable; see Eq. fl62|) . It just happens that in most cases, the irreversible contribution 
turns out to be minimal. However, the reality of the residual entropy is based not on the 
minimal irreversible contribution; it is rather based on the inequality in Eq. ( 160|) . Thus, 
the entropy reduction due to confinement in the unconventional view cannot be attributed 
simply to the presence of irreversibility. Its root-cuse, if it exists, has to be found elsewhere. 

E. Goal and Outline 

Our aim in this work is to follow the recently developed non-equilibrium thermodynamics 



to specifically derive non-equilibrium thermodynamic principles, which s 



to the well-known thermodynamic principles for equilibrium bodies 0, [s], 26|. We must 
ensure that these principles do not violate the second law in anyway. We start by assuming 
the existence and continuity of Sq, and follow I and II. The existence assumption should 
not be surprising as the existence of 5*0 is required to state the second law for an isolated 



lould reduce 
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system Sq; see Eq. ([2]). The equality in Eq. (|2]) occurs only when the isolated system is in 
equilibrium. We will prove here that the existence and continuity of 5*0 implies the existence 
and continuity of the Gibbs free energy [l, 2| 

G(To, Po, t) = E{t) - ToS{t) + PoV{t) (11) 

for the open system S at fixed temperature Tq and pressure Pq of the medium. (Even 
though it is not a common terminology, we will use "open" to denote any body that is not 
isolated in this work.) The continuity of 5*0 is easily traced to the continuity of microstate 
probabilities, as discussed in Sect. HVl The second law for the open system reduces to the 
well-known inequality 

^ < 0; (12) 

the equality occurs when the system achieves equilibrium with the medium; this is exempli- 
fied by the curve ACK in Fig. |2l It is evident from Fig. Hlthat glass free energies Ggl,cv(Tci) 
and Ggl,uv(^o) approach Gsch{To) from above, in conformity with Eq. (fT2|) . Thus, their 
behavior during relaxation is not sufficient to rule out either one as unphysical. On the 
other hand, the continuity of the Gibbs free energy immediately rules out any discontinuous 
entropy reduction due to component confinement as discussed above but does not rule out 
a continuous entropy reduction. 

The continuity of the Gibbs free energy also establishes the continuity of the entropy 
of an open system in internal equilibrium by recognizing that the enthalpy of the system 
remains continuous at the glass transition. The existence of internal equilibrium still allows 
irreversible relaxation in the system. We find that as a consequence of the second law, the 
entropy of the system in internal equilibrium decreases during relaxation, thus validating 
Eq. ( |T3l) and invalidating Eq. ( fT4l) . The latter requires the entropy to increase during 
relaxation as shown by the relaxation arrow for Sgl.uv in Fig- El Thus, the second law 
supports a continuous entropy S during component confinement. On the other hand, the 
Gibbs free energies of the two glasses satisfy the second law in Eq. f|T2|) . In view of the 
incorrect behavior of Sgl,vy, this is a very surprising result. It appears that 5'gl,uv(To) and 
Ggl,uv(^o) do not have the same content. We will find an explanation of this surprising 
statement later. 

The statement about the continuity of the entropy of the system can be made stronger by 
removing the requirement of internal equilibrium by obtaining a statistical formulation for 
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it. However, we only need the above weaker continuity statement to rule out the unconven- 
tional view of the entropy during the glass transition and to support Fig. [2l The statistical 
formulation of the entropy of an isolated system in terms of its microstates is given by Eq. 
(jin]) in compliance with the second law and which satisfies the above thermodynamic princi- 
ples. We then discuss how the entropy is defined for an open system under all conditions in 
terms of microstate probabilities. This derivation also establishes the existence and continu- 
ity of the entropy of an open system by recognizing that the entropy of the medium remains 
continuous during relaxation. This again establishes that continuity of entropy when com- 
ponent confinement occurs. The continuity of the entropy of a body is a generic feature of 
a statistical system, and is formulated as one of the thermodynamic principles. Therefore, 
it should also not come as a surprise that the assumption of the existence of the entropy 
for non-equilibrium states also forms the cornerstone of the conventional non-equilibrium 
thermodynamics jsH]. We then proceed to explain the residual entropy in statistical terms 
and show that there is no reason for it to vanish. We also show that irreversibility alone 
cannot justify any entropy reduction. 

The layout of the paper is as follows. We briefly discuss the two alternative views of 



the glass transition in Sect. HIl which is then 



bllowed by a short review of some of the 
sl in Sect. IIIII Non-equilibrium statistical 



useful concepts and results derived earlier 

m 

entropy is introduced in Sect. |IV] and the concept of residual entropy is discussed. We 
then list and discuss important thermodynamic principles that should be valid for a body 
in internal equilibrium in Sect. |Vl Sect. IVIII is the important section, where we discuss the 
consequences of the non-equilibrium thermodynamics. We discover in Sect. IVII Dl that the 
unconventional view is internally inconsistency. The last section contains our conclusions. 



II. TWO VIEWS OF GLASSES 
A. Conventional Glass View (CV) 

The small amount of irreversible contribution to the entropy of glasses does not mean that 
glasses are in equilibrium (with the surrounding medium); rather, it implies that the fast 
degrees of freedom in glasses have equilibrated, and the slow degrees of freedom, although 
not yet equilibrated, are changing so slowly as to be almost unchanged over the observation 
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time. More details about this point can be found in I. For most glasses, the success of Eq. 
ffTOj) rests on the premise that within 5% of the residual entropy, as noted above, diS ^ 
so that dS ~ dcS. As a consequence, S{0) for a glass in Eq. ( ITU]) is correct within 5%; 
its entire magnitude is only weakly affected by the irreversibility encountered during the 
glass transition. This view will be called the conventional view (CV) of the glass transition. 



Gutzow and Schmelzer 



39| draw an analogy of diS ~ in slowly varying glasses with 



Prigogine's principle of minimum entropy production [ll|]. However, we will not assume 
diS ^ so as to make our argument as general as possible. Later in this work, we will derive 
an almost identical expression to Eq. flTUl) as an identity; see Eq. fIBS]) . In the conventional 
view, the entropy of the glass is schematically represented by the short dashed curve in Fig. 
El the upper solid curve represents the entropy of the supercooled liquid (SCL) that would 
be eventually observed if we wait long enough for the system to equilibrate. The point 
where 5'gl,cv(^o)) branches out of Sscl(Tq) is the glass transition temperature Tog and is 
determined by the choice of observation time robs- As the entropy of the conventional glass 
lies above the supercooled liquid entropy 

'S'gl,cv(^o, ^) > Sscl(To) for Tq < Tog, 



the entropy of the glass actually decreases in time 

dSGL,CY(To, t) 

dt 



< 0. (13) 



The entropy of the conventional glass (blue short dashed curve) approaches > if we 
extrapolate to absolute zero for most glasses. In this view, the glass continuously emerges 
from the supercooled liquid at Tog. Both entropy functions are continuous with only gradual 
variation determined by their respective heat capacities. 



B. Unconventional Glass View (UV) 



Despite the success of the conventional view that provides a justification of the residual 
entropy, there have been several attempts to argue that the residual entropy is not a physical 



concept 



34-: 



It has been suggested that the entropy of a glass should approach at 
absolute zero, regardless of how it is prepared. In other words, the entropy of a glass 
has no memory of its mode of preparation at absolute zero, notwithstanding the fact that 
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FIG. 3: The entropy of the supercooled hquid (5'scl: upper sohd curve) and of the CV-glass 
(S'gl.cv^ short dashed curve) in the conventional view (CV) are shown schematically. We are 
considering an isobaric cooling experiment. The entropy Sgl.cv of the CV-glass continuously 
emerges out of S'scL at Tog and approaches the residual entropy at absolute zero (not shown). 
Schematic form of the entropy S'gl.uv in the unconventional view (UV) is shown by a discontinuous 
entropy loss (thin dashed vertical line of magnitude Sr) at Tog or a continuous entropy loss (dashed 
curve) of 5r from 5scl to the lower solid curve during the glass transition across Tog. Note the 
presence of an inflection point in the dashed piece of S'ql.uv; the latter approaches at absolute 
zero upon extrapolation (not shown). In the presence of an additional variable ^, see text, the 
entropy curves turn into surfaces defined over the Tq — ^ plane. 

most glasses do not represent equilibrium states. The total loss of memory of the mode of 
preparation is the hallmark property of an equilibrium state, while a non-equilibrium state 
remains strongly correlated with its mode of preparation until it fully equilibrates, which is 
reflected in its relaxation towards equilibrium. Despite this, the above suggestion has been 
made repeatedly in the literature. If it is true, it will require either an abrupt or a gradual 
but rapid reduction in the entropy as the supercooled liquid (SCL) turns into a glass (GL) 
at (in the former case) or around (in the latter case) Tog as we cool SCL. The amount by 
which the entropy must decrease during the glass transition (by confinement to one of the 
possible exponentially large number of components) should be comparable to S^. This is 
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the unconventional view (UV) described above. The abrupt entropy reduction during the 
glass transition is shown by the thin vertical dashed line of height at Tog to the lower 
solid entropy curve of the unconventional glass shown as Glass (UV) in Fig. |31 The gradual 
entropy reduction is shown by the dashed red portion with an inflection point, which falls 
very rapidly below the entropy of the supercooled liquid. The total additional entropy 
reduction is equal to S'r. After the drop, the dashed curve connects with the schematic 
lower solid Glass (UV), which when extrapolated to absolute zero will now give a vanishing 
entropy, thus resulting in no residual entropy. We note that because of the entropy reduction, 

SGL,vv{To,t) < S'scl(To) for To < Tog, 

so that the entropy of the glass actually increases in time 

dSGL,uviTQ, t) 
dt 

which is of opposite nature to that in Eq. (fT3|) . 



(14) 



C. Gibbs free energy: VConventional versus Unconventional 

Thus, it appears that there are two competing views at present as far as the entropy of 
the glass is concerned with contradicting behavior of the entropy during relaxation as seen 
from Eqs. f lTSlfT^ . However, they both accept the same entropy function 5'scl(^o) of the 
supercooled liquid before it turn into a glass, and is given by the solid blue curve indicated 
by "Supercooled Liquid." As the enthalpy H and the volume V of the glass are known to 



be continuous during the glass transition 



17H20|], both views agree on the same functions 



Hgl(To) and Vgl(To) for the glass. The Gibbs free energy G(To) has been identified in I as 
G{Tq) = H(Tq) — TqS(Tq) in terms of the entropy of the system. Conversely, the entropy 
can be obtained from the Gibbs free energy. Then the continuity of the entropy in the 
conventional view ensures that the Gibbs free energy Ggl,cy(To) = Hgi^(To) — ToSgl,gy(Tq) 
of the glass is also continuous. This is shown by the short dashed piece in Fig. |4]that emerges 
out of the lower solid curve Gscl{To) without any inflection point. In the unconventional 
view, the Gibbs free energy abruptly jumps upwards by an amount close to T^Sr for an 
abrupt reduction of entropy or rises rapidly along the dashed curve with an inflection point to 
meet continuously with the upper sohd curve to give Ggl,uv(^o) = Hqi^{Tq) — To5'gl,uv(^o) 

16 




FIG. 4: Schematic form of the free energy GscL (lower soHd curve), Ggl.cv (short dashed curve) 
and Ggl.uv (vertical jump at Tog or a continuous dashed piece along with the sohd upper curve). 
The dashed piece has an inflection point A, so that it has a region of positive curvature at higher 
temperatures. There is no inflection point in Gql.cv and GscL- In the presence of an additional 
variable ^, see text, the Gibbs free energy curves turn into surfaces defined over the Tq — plane. 

in Fig. m Both Gibbs free energies coincide at absolute zero: 



G'(O) = Ggl,cv(0) = Ggl,uv(0) = i/GL(0), 



as shown in Fig. |H Otherwise, 



GglMTo > 0) < Ggl,uv(To > 0) for Tq < Tog. 



(15) 



The abrupt entropy reduction in the unconventional view has been recently criticized by 



Goldstein 



4l| . Gutzow and Schmelzer 39|, Nemilov 40|, Gujrat 



ti 0, 



42| . and Johari 



43| . We will not go into their arguments here for which the reader is advised to go to 
the original source; however, Johari 43| provides a good summary of these approaches. In 



addition, Kozliak and Lambert 



44| provide a new look at the issue of the residual entropy. 
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These attempts are different from the approach we take here, although we arrive at a similar 
conclusion. These previous attempts neither consider a continuous entropy reduction nor 
do they visit the issue of how to define the entropy for non-equilibrium states from first 
principles, or how to interpret the residual entropy from a statistical point of view, the 
issues that guide us here. Therefore, the previous attempts do not reveal a pathological 
aspect of the continuous entropy reduction in that the behavior of the entropy 5'gl,uv(^o) 
in Fig. [3] is not consistent with the entropy obtained from the Gibbs free energy Ggl.uvI^o) 
in Fig. m This internal inconsistency poses a major problem for the unconventional view. 

If the entropy can be defined for a non-equilibrium body from first principles and if it 
satisfies the second law, then we should be able to test whether component confinement 
results in an entropy reduction or whether the residual entropy is real for a non-equilibrium 
body. 

III. IMPORTANT CONCEPTS AND USEFUL RESULTS FROM I & II 

Quantities pertaining to Sq are denoted by a subscript 0, for S without any subscript 
and for S by a tilda at the top. The system, although macroscopically large, is considered 
to be a very small part compared to the medium S, so that the system can only create a 
very weak disturbance on the medium. As said earlier, a body refers to any of the three 
systems S, S and Sq. For the sake of convenience, quantities pertaining to a body will be 
denoted without any subscript. 

A. State variables and internal equilibrium 

We first need to clarify the concept of the macroscopic state of a body. As equilibrium and 
non-equilibrium thermodynamics is an experimental science, it must be based on observables, 
which we take to be extensive and collectively denote by X(t), t being the time measured from 
some convenient reference point. Any thermodynamic state function must be a function of 
state variables Z(t) consisting of observable X(t) and internal variable The observables 
are independent continuous variables that can be controlled and manipulated by an observer 
so that their values will allow the observer to differentiate between different macroscopic 
states (macrostates) of the same system. Of the observables, the energy E, volume V and 
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the number of particles N play an important role. The observables are important as they 
remain fixed for an isolated system. The internal variables (see II for details, where they are 
denoted collectively by I(t) = {Ik(t)}) cannot be controlled by the observer, but play a very 
important role during relaxation in non-equilibrium bodies such as glasses. In addition, the 
state function may also be an explicit function of time t. It is useful to replace I{t) by the 
dimensionless vector ^{t) = {Ckit)} defined according to 

_ Ikjt) - 4(oo) 

^'=^'^-4(0)-4(oo)- 

The variables ai'c usually known as "internal order parameters" or the "degree of advance- 
ment." Being a state function, the entropy of a body, whether in equilibrium or not, is a con- 
tinuous function of the state variables and will be expressed as S{Z{t),t) or S(K{t),^{t),t). 

When the body refers to the isolated system Eq, we should remember that its observ- 
ables are fixed so that Xq has no ^-dependence. The entropy of Eq must be expressed as 
S'o(Xo,Io(t), ^) or S'o(Xo,^o(^)' ^) while it is relaxing towards its equilibrium state. Once it 
reaches equilibrium, its entropy can only be described as a function of its observables, so that 
it must be expressed as 5'o,eq(Xo), which is now a constant, independent of t. The explicit 
t-dependence has disappeared. At the same time, Co(^) ^^^^ longer ^-dependent. This 
means that in equilibrium, ^o,eq must become constant, which from its definition above must 
vanish identically: 

^0,eq = 0. (16) 

Moreover, as discussed in II, lo{t) is no longer an independent variable; it becomes a function 
of the observable Xq: 

Io,eq — Io,eq(Xo). 

For an open system in equilibrium, also must become constant in equilibrium; similarly, 
leq is also no longer an independent variable; it becomes a function of the observable Xeqi 

^eq = 0> leq = leq(Xeq). 

The case when the entropy of a body has no explicit i-dependence, so that the entropy 
appears as 

s{x{t),m, (17) 
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is quite special as discussed in I and II. In this case, the system is said to be in internal 
equilibrium. Accordingly, 

under internal equilibrium. The entropy has the maximum possible value for given Z{t). In 
this case, the derivatives 

Y(t) _ dS{Z{tj) A(t) _ dS{Z{t)) 

of the entropy with respect to Z(t) give the field Y(t) and affinity A(t) characterizing the 
body; see II and below. They are defined only when the body is in internal equilibrium or 
in equilibrium. Otherwise, the derivatives do not have any physical significance. 

It will prove useful to think, at least in a formal sense, that the internal variable can also 
be kept fixed. Consider the system in internal equilibrium at the instance when its state 
variables are 

XiE=X(t), IlE=I(t), 

and its field and affinity vectors are 

YiE=Y(i), AiE= A(i). 

The medium E is characterized by Yq, and Aq = 0, and will be expressed as S(Yo, 0); see 
II. We now disconnect this system from the medium and bring it in contact with another 
medium E(Yie, Aje) that is characterized by Yje and Aie. The system is in equilibrium 
with this medium, with the average values of its observables and internal variables given 
by XiE and Ije, respectively. Being in equilibrium, it must satisfy all the requirements of 
stability. In particular, the heat capacity of the system must be non-negative: 

^^y;,,A:, > 0, (20) 

where Y^^ represents all fields except Tie. One can also formally consider this equilibrium 
system as an isolated system with fixed Xje and Ije by disconnecting it from the second 
medium. This isolated system is still in equilibrium, as its entropy is at its maximum for 
given XiE and Ije. The heat capacity of such a system should also be non-negative 

Cx' ,I:e > 0, (21) 
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where Xj^ denotes all observables except energy E{t). If we connect this system to a medium 
which only controls Eis and Vjs of the system by its fixed Tie and Pie, respectively, with the 
remaining observables Xj'p, and Iie kept fixed, then the stability of such a system ensures 
that the following heat capacity is also non-negative: 

Cp T > 0. (22) 
B. Gibbs free energy of the system 

From now on, we will considerably simplify our discussion by considering only E, V and 
as the only observables of a body, and allow only one internal variable in describing it, 
which we denote simply by ^. Moreover, we will keep fixed and allow the possibility 
of fluctuating E and V due to exchange with another body, the medium. Accordingly, we 
will not exhibit A^ in expressing any state function of the body. Thus, X will refer to E, V 
and Z will refer to E, V, ^ from now on. We also assume that the medium is in internal 
equilibrium, so that its entropy S{Z) does not have an explicit t-dependence (so that Eq. 
(ITS]) is satisfied), although it is not in equilibrium with the system. It follows from the 
internal equilibrium of S that the associated fields are 

To 

Here, Zq denotes £"0, Vq and ^0 of ^o- The affinity Aq^ = for the medium. The entropy 
So{t) of Eq is a continuous function of time t, at least twice-differentiable with respect to its 
extensive variables, and satisfies Eq. ([2]). As discussed in I and II, it can be written as the 
sum of the entropies S{t) of the system and S{t) of the medium 

SoiEo, Vo, eo, t) = S{E, V, e, t) + S{E, V, (24) 

by invoking their quasi-independence. We will provide a direct proof of this sum later in 
Sect. IIVI when we identify the entropies as statistical quantities in terms of probabilities of 
microstates. In terms of 

H{t) = E{t) + PoV^(t), G{t) = H{t) - ToSit) = E{t) - To^(t) + PoV{t), (25) 

which are the time-dependent Gibbs free energy and the enthalpy, respectively, of the system 
S at fixed To, Pq and Aq^ = 0, we have 

So{t) -So = S{t) - H{t)/To = -G(t)/To, (26) 
21 
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Zo 



as shown in I. where 5*0 = S{Eq,Vo,^o) is independent of the system. 

Theorem 1 The Gibbs free energy G{Tq, Pq, t) of the system given in Eq. ( fJTl) is determined 
by the temperature and pressure of the medium and is equal to 

GiTo,Po,t)^To[So-Soit)]. (27) 

This is an important observation for the following two reasons. Firstly, it shows that the 
Gibbs free energy is defined even though we may not define the temperature and pressure 
of the system, a situation that occurs when the system is not even in internal equilibrium. 
Of course, we are assuming that the energy, volume and entropy of the system are defined 
under all circumstances. While there is less doubt that this is possible for the mechanical 
quantities (such as energy and volume), the concept of entropy valid under all conditions 
requires some care; see I and Sect. [IV] later. In engineering context, the above Gibbs free 
energy G(To, Po,t) is also known as exergy or availability [45]. 

Secondly, the above conclusion makes an important statement about the continuity of 
G'(t), which we now explain. We first recall that the entropy S'o(t) is a continuous function. 
Now, 5*0 is the entropy of the medium in internal equilibrium (but at Eq,Vq,^q). As said 
above, the properties of a body in internal equilibrium is similar to those of the same body 
in equilibrium, when treated as an isolated system. Since the entropy of an isolated system 
is taken to be continuous, 5*0 must be continuous. This ensures that the difference 5*0 (t) — 5*0 
is also continuous, which then proves that G(To,Po)^) is a continuous function. 

Since the system is an extremely small part of the isolated system, we have 
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Zo 




A, 



0? 



to a high degree of accuracy. Thus, we have 

50(^0, K),6o,t)-5(Eo,\/o,6o) 



Jt L 







dt dt 

where we have used Aq^ = 0, the constancy of Eq, Vq and the second law. The above equation 
is valid even it .^o is not constant, and immediately proves Eq. ( fT2l) . 

Thus, we come to another very important conclusion, which we state as a theorem: 

Theorem 2 The Gibbs free energy G{Tq, Pq, t) of the system is always a continuous function 
of time that continuously decreases in time in accordance with Eq. ( fl^) . 
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The above result does not require any knowledge of the state of the system E. In par- 
ticular, the continuity of G{t) must remain valid even if the processes going on within 
the system are chemical reactions, ergodicity loss, chaos, turbulence, explosion, etc. or 
just normal quasi-static slow processes going on in S. Thus, the continuity of G{t) is 
valid when the system is not in equilibrium with itself or with the medium. As time 
goes on, G(To, P(),t) will continue to decrease towards the equilibrium value G(To,Po) = 
E{To, Po) - ToS{To, Po) + PoV{To, Po). 

We now follow its consequences. Let us assume that the system is in equilibrium with 
the medium at some temperature Ti and pressure Pi of the latter. The Gibbs free energy 
of E in equilibrium is 

G(Ti, Pi) = P(Ti, Pi) - Ti5(Ti, Pi) + PiV{Tu Pi); 

obviously, equilibrium quantities do not have any explicit t-dependence. At t = 0, the 
system is brought in contact with a medium at the current temperature Tq and pressure Pq. 
As the energy, entropy and volume have not had any time to change, the new Gibbs free 
energy at t = will be 

G(To, Po, t = 0) = E{Tu Pi) - To^(Ti, Pi) + PoV{Ti, Pi). 

This is true regardless of what changes may occur at t > 0. It is clear that the two Gibbs 
free energies are very different. Their difference is 

AG{t = 0) = -S{Ti, Pi) AT + V{Ti, Pi)AP ^ 0, (28) 

where Ag = ^current — 9prior for any quantity q. The difference AG(t = 0) at t = 0, however, 
has no physical significance for the second law as the current isolated system with the medium 
at To, Po is different from the previous isolated system with the medium at Ti, Pi, while the 
second law refers to the behavior of a body in time under fixed macroscopic conditions. 
We should also note that AG{t = 0) — > as AT and AP vanish simultaneously for any 
system. Thus, if the macroscopic conditions do no change, the Gibbs free energy undergoes 
no discontinuity at t = 0. 

For t > 0, however, we are observing the same current isolated system, for which Eq. 
(p6|) holds so that the Gibbs free energy G of the system must vary continuously for t > 
regardless of the kinds of processes going on inside So. 
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It is interesting to note that for infinitesimal changes in the temperature and pressure, 
Eq. fl28l) reduces to 

dG{t = 0) = -S{Ti, Pi)dT + V{Ti, Pi)dP, (29) 

a famihar expression from equihbrium thermodynamics. 

If and when the internal equilibrium has been established in S, we can introduce, as 
discussed in I, its instantaneous fields and affinities, such as the temperature T = T{t) and 
pressure P = P{t): 

1 _ dS Pit) _ dS A^it) _ dS_ ^^^^ 



T{t) dE' Tit) dV Tit) dC 
see Eq. ( fT9l) . These are standard relations for the entropy for a body under internal equi- 



librium and are used commonly in non-equilibrium thermodynamics [1, 2, 8-ll|. To ensure 



that the inequality in Eq. (|2]) remains valid during relaxation, i.e. during approach to 



equilibrium, we must have (see also 46| ) 



Tit) ^ To, Pit) ^ Po, A^it) ^ = 0. (31) 

They become identical only when equilibrium has been achieved. Thus, as long as the 
relaxation is going on due to the absence of equilibrium, the two inequalities must hold true; 
see I. Thus, we come to another important conclusions: 

The use of the internal equilibrium is basically in the spirit of Onsager's regression hy- 
pothesis Non-equilibrium relaxation is governed by the same laws as the relaxation of 
spontaneous fluctuations occurring in an equilibrium system. 

We will assume that in an isobaric cooling experiment, 

Pit) = Po, (32) 

which we will refer to as the existence of the mechanical equilibrium for the system. In this 
case, we find that 



dt \T Tq J dt 

It can be shown on general grounds (see I) that in a cooling experiment under mechanical 
equilibrium 
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during relaxation in glasses. This is also seen experimentally. Thus, in such a cooling exper- 
iment 

T > To, (35) 

the equality occurring only when equilibrium has been achieved. Furthermore, in such 
cooling 

dS IdH 

during relaxation at fixed Tq. The relaxation that occurs in the glass originates from its ten- 
dency to come to thermal equilibrium during which its temperature T{t) varies with time; 
recall that we are considering a cooling experiment under mechanical equilibrium. The re- 
laxation process results in the lowering of the corresponding Gibbs free energy, which is a 
consequence of the second law in Eq. ([2]) and is valid even when there is no mechanical equi- 
librium. However, under mechanical equilibrium, the changes in the enthalpy and entropy 
are in the same direction; see Eq. fl36l) . This then gives rise to the following 

Theorem 3 The lowering of G(t) with time in an isobaric experiment under mechanical 
equilibrium results in not only lowering the enthalpy in a cooling experiment, as observed 
experimentally, but also the entropy S{t) during relaxation: 

dS/dt < 0, (37) 

as shown in Fig. O/or the conventional glass (CV) as it approaches towards the supercooled 
liquid during relaxation. 



C. Gibbs fundamental relation 

It follows from Eqs. (jSO]) and ^ that 

dS{t) = ^^dE{t) + P^^dVit) + ^dm- (38) 

This relation is known as the Gibbs fundamental relation. The first law of thermodynamics 
is codified in the following differential: 

dE{t) = T{t)dS{t) - P{t)dV{t) - A^{t)d^{t) (39) 
= TodS{t) - PodV{t) - Asd^s - Ayd^y - A^{t)dat), (40) 
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where the two new and different internal order parameters or degrees of advancement 

_ S{t) - g(oo) _ -y(oo) 
" 5(0) - 5(oo) ' " y(0) - y (oo) ' 

are determined by the instantaneous entropy and volume governed by their corresponding 
"affinities" 

As = -[T(t) - To] [5(0) - 5(cx))] , Av = \P{t) - Po] [V^(O) - y(oo)] , 

thus, ^s, 'Cv pl^-y the role similar to that of internal variables ^, with As, Ay playing the role 
of their respective affinity. The differential of G(t) and -ff(t), see Eq. (125|) . turns out to be 

dG'(t) = -S{t)dTo + V{t)dPo - Asd^s - Ayd^y - A^{t)d^{t). (41) 

and 

dH{t) = TodS{t) + V{t)dPo - Asd^s - Ayd^y - Ag(t)rf^(t). (42) 

For the isobaric case discussed above, we have Ay = 0, but A^ is normally non-zero. 
However, in almost all applications of classical non-equilibrium thermodynamics to glasses 
to date that we are familiar with, even As is taken to be zero 1^, 20]: 

= 0. (43) 

This is equivalent to having 

T{t) = To, (44) 



which we will refer to here as the assumption of thermal equilibrium. From Eq. fl4T|) . we 
see that G cannot vary in time if we have both thermal and mechanical equilibrium and if 
there are no internal order parameters. Thus, there will be no relaxation unless we allow an 
order parameter: 

Theorem 4 Under the assumption of thermal and mechanical equilibrium, internal order 
parameters are required to describe relaxation observed in glasses. 

Except in Sect. IVIH we will not assume mechanical and thermal equilibrium so as to 
remain as general in our discussion as possible. 
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IV. NON-EQUILIBRIUM ENTROPY 



The results in the previous section are based on the thermodynamic concept of the entropy 
for a body, not necessarily in equilibrium. Only the difference in the entropy has any 
meaning as the entropy itself has no unique value, unless supplemented by Nernst postulate 



2l|. T 



lis should be contrasted with the statistical concept of the entropy originally due 
to Gibbs 3| and based on first principles, which provides a unique value of the entropy for a 
given macrostate. Everyone believes that the Gibbs entropy is the correct thermodynamic 
equilibrium entropy j4|. The situation with statistical interpretation of non-equilibrium 
entropy appearing in Eq. (|2]) may not be so clear, which we now analyze in this section. 
The discussion in this section is very general and neither the medium nor the system is 
assumed to be in internal equilibrium. 



A. Fundamental Axiom of Non-equilibrium Thermodynamics 

As thermodynamics is an experimental science, it requires several measurements on the 
body to obtain reliable results. To avoid any influence of the possible changes in the body 
brought about by measurements, we instead prepare a large number A/" of samples or replicas 
under identical macroscopic conditions on which measurements are made. The replicas are 
otherwise independent of each other in that they evolve independently in time. This is 
consistent with the requirement that different measurements should not influence each other. 
The samples are prepared so that the probability of a sample in microstate j is pj{Z(t),t), 
which we simply write as Pj{t); it is a continuous function of the state variable Z{t) and of t 
for a macroscopically large body. We now state the fundamental axiom of thermodynamics 
as proposed in 

Fundamental Axiom The thermodynamic behavior of a system is not the 
behavior of a single sample, but the average behavior of a large number of inde- 
pendent samples, prepared identically under the same macroscopic conditions at 
time t = 0. 

Such an approach is standard in equilibrium statistical mechanics l26|, but it 

must also apply to systems not in equilibrium as it is required for the reliability of mea- 
surements. For non-equilibrium systems, this averaging must be carried out by ensuring 
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that all samples have been prepared with identical history. This is obviously not an issue 
for systems in equilibrium. We refer the reader to a great discussion about the status of 
statistical mechanics and its statistical nature in Sect. 25 by Tolman j26|; see also the last 
paragraph on p. 106 in Jaynes j^S]. This point has been recently reviewed in [s], I23I . 

The average over these samples of some thermodynamic quantity X then determines the 
thermodynamic average quantity X for the body 

w 

Xit)^J2P^{t)X„ (45) 
i=i 

where Xj is the value of X in the jth microstate of the body, and W is the number of its 
distinct microstates. 



B. General Formulation of Entropy: Isolated system 

1. Gibbs Formulation of Nan- equilibrium entropy 

We begin by considering an isolated system Sq, which need not be in equilibrium. Gibbs 



[3I, Tolman [2^ and Rice and Gray [6i] among ot 
entropy using the Gibbs formulation; see also 



lers discuss at length the non-equilibrium 



231]. It is given by the negative average of 



7] = \np, 

n 

what Gibbs [3] calls the index of probability: 



So{t) = -\np = -ZPa{t) lnp„(t), EPa(t) = 1, (46) 

a a 

where Pa{t) is the probability of the ath microstate at time t. It is shown below that 
it follows from Eq. (145|) . As the microstate probability Pa{t) is a continuous function of 
its arguments Zq and t, the entropy So(t) = S'o(Zo,t) is also a continuous function of its 
arguments Zq and t. It is straightforward to establish [sl, that this entropy satisfies 
([2]). The identification of the entropy with the negative of the Boltzmann if-function [ 
see p. 561], the latter describing a non-equilibrium state, should leave no doubt in anyone's 
mind that the Gibbs formulation of the entropy can be applied equally well to an equilibrium 
or a non-equilibrium system in isolation. 



Lq. 



26|, 



28 



For a generic body, we will write the continuous entropy as 

S(t) ^ -h^ ^ -J^P^it) \np,it), ZpAt) = 1, (47) 
j ^ 

where Pj{t) is the probability of the jth microstate of the body at time t. The justification 
for the entropy expression in Eq. p7|) for an open system, its continuity and its connection 
with the second law will be given below. 

There are two different ways to understand the above entropy formulation. 



2. Ensemble Interpretation- Dynamics being Irrelevant 

We now prove that the entropy is a statistical average given in Eq. fj45|) . We consider a 
large number Af = CWqCKq) of independent replicas or samples of the isolated system, with 
C some large constant integer. Let there be N'a samples in the microstate a. The probability 
of a sample in microstate a is then 

P„ = ATJAf. (48) 
The ensemble average of Zq over these samples is given by 

_ I l^o(Xo) VKo(Xo) 

Z0 = T7 E A/'aZoa= E Pa'Z'Oa, (49) 

a=l a=l 

where Zoa is the value of Zq in the microstate a. This is identical to the statistical average 
given in Eq. (H5l) . 

The number of ways W to arrange the A/" samples into Wq distinct microstates so that 
there are JVa samples in microstate a is given by 

a 

Taking its natural log to obtain an additive quantity 

S = \nW, (50) 

and using Stirling's approximation for the factorials, we see easily that S per sample can be 
written as 

S Wo{Xo) 

77 = - E Pa In Pa, 
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where pa is given in Eq. (HHj) . Thus, S/M is the ensemble average of —Inpa, thus, it is 
nothing but the entropy 5*0 of the body given in Eq. fH6l) . Again, the continuity of S/J\f 
follows directly from the continuity of Pa- 

The maximum possible value of the entropy for fixed Xq is 

S-cmaxlXo) =lnl^o(Xo), (51) 

which occurs if and only if all microstates are equally probable: 

p^{t) ^ l/W^o(Xo), Va G To. 

This maximum value is the Boltzmann entropy of the system. This certainly occurs as 
t — 7- oo for most systems. Then the system is said to be in equilibrium. It is evident that 
the Gibbs formulation in Eq. fH6|) is more general than the Boltzmann formulation in Eq. 
f lSTj) . as the former contains the latter as a special limit. In equilibrium, Af^ = C for all a, 
so that Pa = I/Wq. 



3. Temporal Interpretation- Dynamics being Relevant 

Another way to interpret Eq. fj46l) is as follows, which is quite standard, at least in de- 
veloping the kinetic theory of gases. We consider the time-evolution of an isolated system 
[sl, [3], starting at t = from some initial microstates ao. Then Pa{t) at some later time 
t represents the frequency with which the ath microstate has occurred during this time 
interval. This is the temporal definition of the probability 4|, the traditional way of intro- 
ducing this probability. All these microstates belong to the slice Tq of the microstate space, 
consisting of fixed Xq, but different .^o(^)- At t = 0, Pa(0) = 5a,ao , where 6 is the Kronecker 
delta, and 5*0(0) = 0. With time, this entropy will continue to increase until it reaches its 
maximum for fixed Xq. Let lyo(Xo) denote the number of distinct microstates in the slice 
To. 

For whatever reasons, if it happens that the system is confined to, or has visited during 
certain time interval r (such as the observation time Tobs), only a part Tq C Tq of the above 
microstate space slice, then 

Pair) = 0, Va ^ T'q, (52) 

even though these microstates corresponds to the same Xq. This happens because these 
microstates have not been visited yet. Let Wq^Xqjt) < Wq^Kq) denote the number of 
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distinct microstates in Fq. Then, the conventional entropy in Eq. (H6|) with the above 
condition in Eq. fl52l) would be strictly bounded 

S'o(Xo,r) = - E Pa{t)\np^{t) < lniy^(Xo,r); 

the equality occurs if and only if all these microstates in Fq happen to be equally probable: 

p«(r) = l/W^^(Xo,r), WaeT',. (53) 



Again we notice the generality of the Gibbs formulation in Eq. P6j) over the Boltzmann 
formulation in Eq. (!5T|) . 

The ensemble approach is also very important from an experimental point of view. At 
high temperatures, where dynamics is very fast, it is well known that it agrees with the 
temporal formulation. However, at low temperatures, where dynamics becomes sluggish as 
in a glass, the temporal entropy can be misleading. However, the ensemble entropy still 
gives the correct value. This point has been discussed in 8|. 

C. General Formulation of Entropy: Open system 

Let us now consider our system S with at least one extensive observable that is fixed 



22], which we take to be the number of particles A^. As usual, the system is a small part of 
the isolated system Sq. We wish to determine the entropy of S in terms of its microstates, 
which are indexed by i. Theses microstates correspond to all possible allowed energies and 
volumes. We use 5 to denote the microstates of S with a fixed number of particles A^. A 
specification of the microstates i and a gives a unique microstate specification a. Hence, the 
number of microstates Wq of the Eq is the product WW, where W and W are respectively 
the number of microstates of S and S. Because of the smallness of S relative to Sq, which 
results in the quasi-independence of the system and the medium to a very high degree of 
accuracy, we have in terms of their probabilities 

As usual, these probabilities are continuous functions of their arguments. Now, using 
\n.pi{t)pa{t) = In pi{t) + \iapa{t), and the sum rule 

a i 
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we find that 

i a 

where the two terms in the above equations represent the entropies of the system and the 
medium 

S{t) = -EP^{t) Inp.(t), S{t) = -EPait) lnj95(t). (54) 

i a 

Note that our entire discussion never puts any restriction on the kind of actual processes 
going on during approach to equihbrium. Hence, whether there is any glass transition, any 
loss of ergodicity loss, etc., has no bearing on the formulation of entropies in Eq. (154|) or 
our other general results presented in the previous section. 

Let us imagine isolating the medium by disconnecting it from the system. Its entropy 
is given by S(t) = S{Z(t),t). The isolated medium is similar to the isolated system Eq, 
whose entropy was assumed to be continuous; see the discussion immediately after Eq. (|2]). 
Thus, S(t) must also be a continuous function. The continuity of So(t) and of S(t) then 
immediately leads to the continuity of the entropy S{t) of the system. This result is also 
consistent with the continuity of microstate probabilities Pi(t). The continuity of S(t) is 
independent of the kinds of processes that are going on within it. This discussion then leads 
us to the following important theorem: 

Theorem 5 The entropy S{Z{t),t) of an open system is a continuous function of its ar- 
guments under all conditions. In particular, it is also continuous when the medium is in 
internal equilibrium. 

The above derivation not only justifies Eq. ([T]) or Eq. (12^ for an isolated system, but 
also identifies the proper non-equilibrium entropy for our open system S, which is always 
a continuous function of its extensive arguments and t. The expression for S in terms of 
microstate probabilities in Eq. (15^ is identical in form to what one would use for an 
equilibrium system [J] , except that the microstate probabilities are now explicit functions of 
time t. This entropy then determines the Gibbs free energy and enthalpy in Eq. (!25l) of the 
system. The second law in Eq. ([2]) for the isolated system Sq is now expressed in terms of 
the Gibbs free energy for the system S and is given in Eq. ( fT2l) . provided that the medium is 
taken to be in internal equilibrium. The fact that E is in a medium E is reflected in the fact 
that the G{t) and H{t) depend on To,Pq of the medium, and not on T{t),P{t) and A^{t) 
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of the system. Recall that the system may have no well-defined temperature, pressure, etc. 
unless the system is in internal equilibrium. This is important to emphasize as the validity 
of the second law cannot depend on the establishment of internal equilibrium in the system. 
Thus, the entropy formulation in Eq. (jSlD is very general. 

Notice that the ensemble interpretation of the entropy in a body is due to the choice 
of samples with the probability pf, the dynamics within the body never appears in the 
discussion. Thus, such an interpretation is quite useful when the dynamics in a system 
becomes very sluggish, such as in a glass. 



V. PRINCIPLES OF NON-EQUILIBRIUM THERMODYNAMICS 

Based on the discussion above and the Fundamental Axiom, we can formulate four 
fu„da,„e.ta. pnnc.p.es of „o„-e,un,b.u,„ thennodynaoics 0. These p.ncip.es a.e we.. 

known and accepted in equilibrium thermodynamics. We now discuss why they should also 
be applicable to slowly evolving metastable states such as glasses in which we are interested 
in this work. We will assume that the system of interest is in internal equilibrium, but not 
in equilibrium with the medium. 

1. The first one is the principle of additivity, according to which the total entropy or 
other extensive quantities can be obtained by a sum over different macroscopic parts 
of the body. Each part must be large enough so that the usual argument that their 
surface effects can be neglected as thermodynamically unimportant is valid so that the 



parts become quasi-independent 



This principle is consistent with what one must do for any measurement, which requires 
verification by performing the measurement many times over on different samples and per- 
forming an average over these samples. Each measurement must be performed under iden- 
tical macroscopic conditions. We now argue that different quasi-independent parts of the 
body represent different samples of the body prepared under identical conditions. Consider 
some body S' with A^' particles, and imagine a much larger body S that can be divided 
into a large number of parts S'^, A; = 1, 2, ■ ■ ■ , each part with A^' particles so that they are 
the same size as the body S'. Let us imagine S itself to be a very small part of an isolated 
system Eg, so that various S'^ and S do not affect the internal equilibrium of the much 
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larger medium S. In particular, Tq,Pq and Aq^ = of S are unaffected by and S. We 
can imagine each to be in a medium S^, where the latter is composed of the medium E 
and the remainder of the system S obtained by taking out the part S'^ under investigation. 
Because of the smallness of S relative to S, the new medium also has the same Tq, Pq and 
Aq^ = as the medium S. Thus, each part S'^ experiences the same Tq,Pq and Aq^ = 
as any other part: all parts experience the same macroscopic conditions. Because of this, 
different parts S'^ are nothing but the preparation of the body S' many times over under 
identical conditions specified by the same Tq, Pq and Aq^ = 0. Thus, averaging any quantity 
over different samples is equivalent to averaging over different parts J^'f, of S. This principle 
of additivity must apply to any body even if it is not in equilibrium. Consequently, the 
value of any extensive quantity over the entire body must be a sum over its various parts, 
regardless of whether the body is in equilibrium or not. 

2. The second principle is the principle of reproducibility, according to which the en- 
semble average is equal to the average of the experimental values, also called the 
thermodynamic average. 

This principle follows from the FUNDAMENTAL Axiom according to which thermodynam- 
ics, whether equilibrium or time-dependent, requires several measurements on the system to 
obtain reliable results. To avoid any influence of the possible changes in the system brought 
about by measurements, we can instead prepare a large number of samples under identical 
macroscopic conditions. This is consistent with the requirement that different measurements 
should not influence each other. Then the principle follows immediately from Eq. ( H9|l . 

3. The third principle is the principle of continuity and uniqueness , which states that 
the Gibbs free energy is continuous, and hence a unique (single-valued) function of its 
arguments. 

For equilibrium states, this principle is certainly valid. We extend it non-equilibrium 
systems. The principle follows from Conclusion 2. As said there, its validity is independent 
of the kinds of processes going on within the system S. In the absence of uniqueness, 
different experimentalists will have no way to effectively communicate their results and no 
scientific investigation can be carried out. This is a very important requirement if the second 
law for an open system in the form of Eq. ([T2|) has to make any sense. The existence of the 
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Gibbs free energy, as demonstrated in Sect. IIIIBI does not require the system be at least in 
internal equilibrium; it only requires the medium to be in internal equilibrium. Therefore, 
it is most certainly also valid when the system is in internal equilibrium. 

4. The last principle is the principle of stability, according to which the heat capacity, 
compressibility, etc. must remain non-negative for the system to remain stable. 

We have already discussed the non-negativity of the heat capacity in Eqs. (I20tl22p . Similar 
conditions also apply to compressibility and other response functions. 



VI. ENTROPY AS A MACROSTATE PROPERTY AND COMPONENT CON- 
FINEMENT 

Can there anything be wrong with the original argument that suggests entropy reduction 
due to confinement as noted in Sect. ^ After all, if the system is confined to just one 
component, should not the confinement entropy be zero? What can be wrong with such 
a simple argument? To answer this question, we need to understand the concept of tjbe 
statistical entropy, to which we now turn. More details can be found in a recent review 



A. Non-interacting Ising Spins: Ensemble and Temporal Descriptions 

Let us consider an ensemble of non-interacting Ising spins located at the sites of a 
lattice with sites. All possible W = 2^ microstates have the same energy E = 0, and 
each microstate Ij has the same a priori probability pj = 1/W. The spin macrostate I, on 
the other hand, is only specified by the number of up and down spins, and not the sequence 
of spin states on the lattice. One can relate this problem to picking balls with replacement 
from an urn containing equal number of balls of two colors. After picking a ball, we must 
place it back in the urn. The sequence of the colors of the balls in time with replacement 
determines a microstate Bj of the A^ balls. A ball macrostate B is determined by specifying 
only the number of balls of each color, but not the actual sequence of colors. We can identify 
the state of the kth spin in Ij with the color of the ball picked at the kth attempt in Bj. 
Then the two problems are identical, except that we are considering an ensemble of the Ising 
spins, whereas we are considering a temporal description for the balls. Indeed, in the latter 
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case, we have a temporal evolution of the state of a single ball in time. Now, there is no 
dynamics that allows a ball microstate Bj to change into another ball microstate Bj. There 
may be some dynamics that could change a spin microstate Ij to another microstate Xj. 
Whether there is any microscopic dynamics specified or not is irrelevant in determining the 
entropy for the spins. 

The entropy for both systems, spins and balls, is 

S = N\n2. 

It is also the maximum entropy, indicating that we are dealing with an equilibrium state. 
This entropy will be the same even if there were no dynamics changing the spins, such as at 
absolute zero in classical mechanics. When a system is confined to one of disjoint components 
of the phase space from which it cannot escape, we also encounter a situation with the 
absence of a djTiamics, the dynamics which takes the system from one component to another 
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23| . Accordingly, every realization of the balls or spins will remain in its microstate 
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m 



forever, just like a glassy system which remains confined to a disjoint component 
0,0] due to kinetic freezing. However, as the Fundamental Axiom states, the entropy is 
an average quantity obtained by an average over all samples or microstates, as shown in Eq. 
(H7|) . This is equivalent to saying that the entropy is determined by the macrostate, which 
represents a collection of microstates, each with certain a priori probability. The entropy 
has a contribution from all of these microstates. It is not the property of a single microstate. 
This point should be very clear form the derivation of relating S/N' in Eq. (!50|) with the 
entropy. The evaluation of S is due to the choice of various samples; the dynamics is not 
part of the derivation. As a consequence, the entropy is unaffected by whether microstates 
have any dynamics for change or not in time. One only needs to know the probability 
distribution for the microstates. 

For a kinetically frozen glass, different samples will correspond to a glass confined to any 



one of the components 
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23j. We have no information as to which component a sample 
is frozen into. Thus, the entropy must be obtained by averaging over these different samples 
or disjoint components. We cannot just consider one particular sample, since the entropy is 
not the property of a single component. 
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B. Probability Collapse and Entropy Reduction 

Is it possible to justify that the entropy is determined by considering just one of the 
components in which the system is confined? If such a determination is possible, it must be 
because the entropy discontinuously decreases when the component confinement occurs. 

The reduction in entropy is most certainly possible if we perform a measurement to iden- 
tify the microstate in which a given sample is. To identify a particular microstate requires 
complete information about the body, which is ordinarily unfeasible. In order to identify 
any particular microstate, we need to perform a very special kind of measurement, which 
we will call a microstate measurement [8^], that provides us with the complete information 
about the body in its current microstate jo- For the Ising model, this requires determining 
the orientations of each of the spins. After the microstate measurement, we know with 
certainty which microstate the system is in. Before the measurement, a given sample is 
known to be in one of the microstates. The probability po that the sample is in microstate 
jo is Po = l/W. This probability changes discontinuously from po = 1/W to po = 1 imme- 
diately after the measurement. The effect of the microstate measurement is to also reduce 
the probabilities of all other microstates j' 7^ jo to pjf = for this sample. Thus, 

Pj ~ ^j,jo^ 

immediately after the measurement. We will speak of probability collapse to indicate this 
change in the probability brought about by the microstate measurement in this work. The 
entropy also vanishes in an abrupt fashion immediately after the measurement from the 
initial value of In W in accordance with complete certainty about the body. 

While quite an appealing argument for the justification, it overlooks two important facts: 

1. In experimental glass transition, no such measurement is ever made that identifies 
precisely which component the glass is frozen in. Such a measurement will tell us 
precisely the positions and momenta of all the particles which then allows us to 
decipher which particular component the glass is in. 

2. Because of the lack of such a measurement, we must determine the entropy by averag- 
ing over all components as shown in Eq. (!55|) by considering all possible samples with 
the probability distribution px for the components. 
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C. System Confined in a Component: The Residual Entropy 



In real glass formation, no probability collapse occurs as no microstate measurement is 
ever performed to identify the particular component the sample is in. The sample still has 
the same probability to be in any of the C components it had before the confinement occurs. 
The confinement does not alter the probability distribution. All we know for sure is that a 
glass sample is in any one of the C components. We do not know the actual component it is 
in. This situation is identical to the die which is known to be in one of the six possibilities 
when it is concealed by the cup, or to the non-interacting Ising spins discussed earlier. As 
long as the die is concealed, the probability of an outcome of an unloaded die remains 1/6, 
regardless of what is on the top face. Even if the outcome cannot change when we lift the 



Is 



for more 



cup, the outcome is most certainly not certain when concealed by the cup; see 
details on this point. The information that the glass has been formed is most certainly 
not equivalent to knowing precisely the particular component in which the glass is trapped. 
The latter will require the collapse of the original probability distribution by a microstate 
measurement. The entropy is obtained by taking the average over all the samples, and not 
only one glass sample, in accordance with the Fundamental Axiom. 

Let us now follow this line of reasoning. The entropy S{t) can be written as a sum over 

all the components, indexed byA = l,2.---,C 

c c 

A=l jx a=l ja 

where jx represents one of the microstates in the component A, the first sum is over all the 
components and the second sum is over all microstates in each component. The above way 
of writing is an identity and applies under all cases. Thus, it applies whether confinement 
occurs or not. Introducing 

c 

Pa = J2p^^^^P^ = ^' 

jx A=l 

we can rewrite S{t) as follows: 

c 

Sit) = J2PxSxit) + Scit), (55) 

A=l 

where 

^ Px P\ 

Jx 
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denotes the component entropy in a given component A, and 

c 



denotes the confinement entropy due to various components. For an unbiased samphng of 
the components, Sc{t) must be at its maximum, which requires 

Px = l/C. (56) 

For the case of the die, = ln6 and Sx(t) = to give S' = ln6, as noted earher. 

In real glasses, the entropy S\(t) vanishes as the glass approaches absolute zero. In this 
case, the entropy at absolute zero reduces to Sc{t), which is what is customarily called the 



residual entropy S'r [19|, l20(| . We thus define the residual entropy as 

c 

Sk = -"^Pxlnpx; (57) 

A=l 

when the system explores various components without any bias so that Eq. ( !56|) holds, the 
residual entropy reduces to 

Sn = \nC. (58) 



D. Residual Entropy of Subsystems and Calorimetric Measurements 

Another way to understand Eq. ( l57l) is to recall the additive of the entropy. We imagine 
dividing the system S into several macroscopically large but equal parts of size , these 
parts representing many samples of a smaller system S'. Each part now represents a glass 
trapped in a component A' of the smaller system S'. Let C denote the number of disjoint 
components for E'. Then the number of disjoint components C for E is given by 

c = {cY/^' , 

where N/N' denotes the number of S' parts in S. Again, we cannot be sure of which 
component A' each part is trapped in. The components appear with probabilities py for any 
part so that 

C 

A' = l 
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for each part. The probabihty of S is obtained by adding this entropy over all parts so that 

If the components appear with no bias, then 

'^R=^lnC', (59) 

which coincides with the result in Eq. ( 158|) for S. The additivity principle clearly shows 
that the residual entropy for S is not zero. 

It is evident now that to conclude that the residual entropy has vanished just because a 
sample has frozen into a single basin or glass form is incorrect. The entropy reduction for S 
only happens if a microstate measurement is performed to identify the particular component 
the E-glass is in. This discussion also shows that calorimetric measurements explore different 
glass components associated with the subsystems, so that they reveal a non-zero residual 
entropy. 



E. Spontaneous Symmetry Breaking versus Glass Confinement 

It is relevant at this point to contrast the above component confinement in glasses with the 
idea of confinement that occurs in spontaneous symmetry breaking such as a ferromagnet. 
In the latter case, there exists a symmetry breaking field, whose presence picks out one of 
the components, for example A = Aq, for the system. This is equivalent to our notion of the 
microstate measurement, except that in this case it is really a "component measurement" 
resulting in picking the component Aq: the application of the symmetry breaking field forces 
the system to be in a particular component Aq. In this process, the entropy of the system will 
be reduced by Sc (note that we have suppressed the time- dependence as one usually considers 
equilibrium situations in symmetry-breaking) due to the probability collapse discussed above 
on p. [371 all thermodynamic averages are determined by the particular component 
Aq. It usually happens that the value of Sc in spontaneous symmetry breaking is not an 
extensive quantity, so the effects of the confinement on the entropy become irrelevant for 
a macroscopic system. There is no entropy reduction per unit volume in the process of 
confinement. However, the effects of the confinement on some thermodynamic quantities, 
such as the order parameter associated with the symmetry breaking, become very important. 
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For glasses, the residual entropy Sc(t) is found to be an extensive quantity. Therefore, it 
is very important to know if there could be some entropy reduction due to confinement in 
glasses. To this date, no one has identified any physical symmetry breaking field analog for 
glasses. Thus, there does not seem to be a way to prepare a glass in a particular component. 
When a glass is prepared, we have no way of knowing which component it is trapped in. 
Hence, one must consider all the components in obtaining any thermodynamic average, such 
as the entropy as we have done above. In particular, there is no entropy reduction per unit 
volume in a glass transition. 

We have discussed the unsuitability of time-average (instead of the ensemble average) 
elsewhere Q|, and we refer the reader to this for more details. We will only make the following 
brief comment here. It happens that at low temperatures the approach to equilibrium takes 
more time than feasible due to experimental constrained. Most measurements last a short 
period of time. The temporal average over an extended time period has nothing to do with 
information obtained in measurements that may take a fraction of a second or so. Unless 
the system is already in equilibrium at the time of the measurement, different measurements 
carried out on the system at different instances will give different results; the system has 
a memory effect, when the system is not in equilibrium. Thus, temporal average is not 
desirable for glasses. 

In contrast, the ensemble average provides an instantaneous average and thus bypasses 
the above objection of the finite measurement time. 



F. Role of Irreversibility on the Residual Entropy 

There is another way to understand the inequalities in Eq. ([6]). In a vitrification experi- 
ment from a state A at temperature Tq in the supercooled liquid state which is still higher 
than the glass transition temperature to the state Aq at absolute zero, we have along the 
path A-)-Ao 

^(0) = S{n) + Jd,S + JdiS, (60) 

A A 



where we have set A5'c = as there is no latent heat in the vitrification process, and where 




dS = deS + diS pj, l2|, l9l-|ll|. with diS > representing the irreversible entropy generation 



and dcS = CpdTo/TQ; see Eq. ([3]). Since the second integral in the above equation is always 
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non-negative, we obtain 

3(0) > Se.pt(0) = Sin) + JCpdTo/To, (61) 

To 

in accordance with Eq. as expected. The forward inequahty is due to the irreversible 
entropy generation. Thus, the entropy at absolute zero must be larger than or equal to 
the right hand side 5'expt(0), which can be determined by performing a cooling experiment. 
We take Tq to be the melting temperature Tm, and uniquely determine the entropy of the 
supercooled liquid at Tqm by adding the entropy of melting to the crystal entropy 5'cr(Tom) 
at Tqm- The latter is obtained in a unique manner by integration along a reversible path 
from To = to To = Tom: 

Scr{Tm) = Scr{0) + J Cp,cRdTo/To, 



here, S'cr(O) is the entropy of the crystal at absolute zero, which is traditionally taken to be 
zero, and Cp^cr(To) is the heat capacity of the crystal. This then uniquely determines the 
entropy of the liquid to be used in the right hand side in Eq. fl61l) . If experiments show that 
5'expt(0) greater than zero (or greater than S'cr(O) if the latter is not taken to be zero in 
accordance with Nernst postulate), then the entropy S{0) at absolute zero on the left side 
itself must be greater than zero (or greater than S'cr(O)). We will assume that 5'cr(0) = 0. 
Thus, the experimental determination of the right hand side of Eq. ( I6TI) is sufficient to 
unequivocally determine whether S{0) > S'cxpt(O). The inequality in Eq. ( |6T1) takes into 
account any irreversibility during vitrification. As is well known, S'cxpt(O) is found to be 
non- negative, as discussed in Sect. II A[ Thus, there cannot be any doubt that 

Sr = S{0) > 5expt(0); (62) 

experiments invariably give a non-zero value of S'expt(O), which then proves immediately that 
in thses cases, the residual entropy cannot vanish. 

VII. NON-EQUILIBRIUM ENTROPY AND GIBBS FREE ENERGY 
A. Consequences of the Second Law 

The important conclusions from the previous discussion can now be summarized as fol- 
lows. The conclusions are derived under the assumption that the entropy 5*0 of the isolated 
system is a continuous function of Eq,Vo,No, ^o{t) and t. 
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(A) Medium under internal equilibrium 



1. The Gibbs free energy G of the system is continuous function of its arguments 
and decreases continuously during relaxation toward its equilibrium value. 

2. The entropy S of the system is continuous function of E{t), V{t),^{t) and t, which 
follows from Theorem [5l 

(B) Medium and system under internal equilibrium 

1. The entropy 5* of the system is continuous function of E{t), V{t), and ^(t). This 
is a weaker form of the continuity condition, and follows from the continuity of 
the Gibbs free energy in Theorem |2J To see that we write 

S={H-G)/T,, (63) 

and use the experimentally observed continuity of H during the glass transition. 

2. The entropy continuously decreases during relaxation toward its equilibrium value 
according to Eq. fl57|) . 

It follows from either (A2) or (Bl) that the entropy cannot undergo any discontinuity at 
a glass transition. Since the continuity of G follows from the second law, any discontinuous 
entropy reduction will result in the violation of the second law, a point already made by 



Goldstein 



4l| 



It follows from (B2) that the behavior of entropy during relaxation only supports Eq. 
f fT3|) and not Eq. f|T4|) . In other words, only the conventional view of the entropy during 
component confinement can be substantiated by the second law. However, this does not 
prove that the Gibbs free energy G'gl,uv(^o) is incorrect as its variation in time conforms 
to the expected variation in Eq. f lT2|) . The Gibbs free energy Ggl.cvI^o) also conforms to 
the expected variation in Eq. ( !T2|) . Does it mean that Ggl.uvI^o) is just as acceptable as 
Cql.cvI^o) from the point of view of the second law and non-equilibrium thermodynamics, 
even if 5'gl,uv(^o) is not? Or is it possible that both Gibbs free energies are unacceptable? 
Then what should be the correct form of the Gibbs free energy for a glass? Is it possible 
that Ggl.uvITo) is acceptable, even though 5'gl,uv(^o) is found to violate the second law? 
Will this lead to some thermodynamic inconsistent in some way? Is it possible that only 
Ggl,cv(2^o) is acceptable? We turn to this issue now. 
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B. Non-equilibrium Entropy evaluation 



The Gibbs free energy Ggl.cvITo) is shown by the short dashed blue curve which grad- 
ually connects without any inflection point to Gscl(Tq), while Ggl,uv(^o) is shown by the 
solid red curve and by its discontinuous or rapid fall to Gscl(^o) in Fig- IH the latter pos- 
sesses an inflection point A in the dashed red portion. We rule out the discontinuous jump in 
Ggl.uvI^o) as physically unacceptable because of (Al). Thus, we only consider Gql.uvI^o) 
with a continuous fall to Gscl(Tq). 

All known analyses of experimental data have been carried out under the assumption of 
mechanical and thermal equilibrium. Therefore, we will also make this common assumption 
in this section, according to which 

As = Ay = 0. (64) 

Accordingly, we will assume that the instantaneous temperature and pressure T{t),P{t) of 
the system are no different from the constant temperature and pressure Tq, Pq of the medium. 
As there is no longer any need to distinguish the two, we will use the common notation T, P 
for both, with the implicit assumption that they have no longer any t-dependence. 

In our discussion below, we do not specifically endorse the statistical formulation of the 
entropy in Eq. f l5^ . We simply use the thermodynamic notion of the entropy. The resulting 



differential of the Gibbs free energy is now given by 19|, |20| 



dG = -SdT + VdP - A^dC (65) 

The irreversible entropy generation is captured by the presence of ^ in the above equation. 
This form of dG is exactly what is postulated in the traditional non-equilibrium thermo- 



dynamics |9Hl 1| . where the existence of non-equilibrium entropy is postulated. How it is 
actually calculated is irrelevant for the general argument below. This point should not be 
forgotten in the following discussion as the general conclusions should be confirmed by any 
formulation of the entropy. 
The entropy is given by 

S = -{dG/dT)p^^, (66) 

which generalizes the conventional definition of the entropy for an equilibrium system to a 
non-equilibrium system. The above relation shows that we should be able to extract the 
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entropy S from the knowledge of the Gibbs free energy G. Now, it is easy to show [l9| that 
the heat capacity is given by 

Cp,^(T) = T {dS/dT)p^^ = -T{d^G/dT^)p,i- > 0; (67) 

this is in accordance with Eq. f l22|) . We can now integrate Eq. fl67|) to obtain the entropy 
in terms of the heat capacity Cp^^ : 

S{T, = SiO, + + A5c, (68) 



where 5'(0,^) is the residual entropy at absolute zero, and where 

represents the amount of heat at constant P, S,- The Eq. fIBS]) should be compared with 
Eq. ( TTU]) . which is usually criticized as unreliable due to the glass transition. In contrast, 
Eq. (1681) is an identity in our non-equilibrium thermodynamics jl| under the assumption 
in Eq. flM|) . Any irreversibility in the system is now captured by the time dependence of 
^. However, the actual form of the variation is not relevant for our discussion in this work. 
The integration in this identity must be carried out along constant ^. This heat capacity 
cannot vanish at any positive temperature. It can vanish at absolute zero without creating 
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22l |. Any state with 



any conceptual problem. We have discussed these issues elsewhere 
a negative heat capacity will be identified here as unphysical. 

C. Consequences of ^ for Glasses 

Let us try to understand the consequence of the additional variable ^. In the presence 
of ^, all the the curves including the broken red curve in Figs. [3] and H] turn into surfaces 
defined over the T— ^ plane. In particular, the inflection point A that is present in the dashed 
piece in this figure turns into a line A(^) of inflection points. Thus, if we consider a slice 
of these surfaces at a constant C, = ^c, then this slice will give three curves similar to those 
shown in Figs. |3]and|ll with S'gl.uvI^o) and Ggl.uvI^o) still containing inflection points. In 
particular, we will find a region curving down at higher temperatures and a region curving 
up at lower temperatures in Fig. [3l just as we see in the figure for 5'gl,uv(^o)- Nevertheless, 
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the form of the entropy 5'gl,cv(To) or S'gl,uv(^o) for fixed ^ = is monotonic as shown in 
Fig. 131 so that 

^^>0 (69) 



OT J 



for both of them. Since we are allowed to take any value in the above discussion, there 
is no harm in showing the fixed value as a general ^ in the above equation. The above 
inequahty corresponds to Cp^^ > 0, in accordance with Eq. (j67l) . Again, we show ^ and not 
.^c for the reasons alluded to above. In an isobaric experiment at fixed P, the curves shown 
in Fig. [3] need not correspond to a fixed ^ in the glass transition region. In general, we 
expect ^ to vary with the temperature, pressure and the history of cooling. But the point 
to stress is that the presence of the inflection point A in Fig. [3] has no effect on the sign of 
the slope {dS/dT)p^. It remains non-negative for both entropies 5'gl,cv(^o) or 5'gl,uv(^o)- 



D. Forms of Gibbs Free Energies 

The enthalpy Hqi^{T,C,) is the same in both views 

because of the continuity of the enthalpy during the glass transition. The Gibbs free energy 

Ggl{T, = Hgl{T, - TSoLiT, 
thus, has the same value Hqi^{0,C,) in CV and UV at absolute zero: 

G'gl(0,O = G'gl,cv(0,O = G'gl,uv(0,O = ^gl(0,O, 

regardless of the value >S'gl,cv(0,O 'S'gl,uv(0,O take at absolute zero (we obviously do 
not consider the unphysical case of an infinitely large S{T,C,) for a finite but macroscopically 
large system) , depending on whether we follow C V or U V. This is the same as Eq. ( 1T5|) . The 
SCL Gibbs free energy G'slc(^) above Tg = Tog is the same in both views. If it is possible for 
the "Supercooled Liquid" curve to be extrapolated to absolute zero, as shown by the solid 
blue curve in Fig. HJ then Gscl{T) would have the value GsclIO) = ifscL(O) < -f^GL(0,^) 
there. Observe that 

GsLciT) < GGL,cy{T,^) < Ggl.uvI^, 
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for T < Tg. The value Gscl(O) should be strictly lower than GGh{0,$,), as Ggl,cv(^, 
GGh,vv{T,^) is expected to reach Gscl(T) from above during relaxation. The curve Gslc(^) 
presumably approaches absolute zero with a slope equal to zero (5'scl(0) = 0). The entire 
supercooled liquid curve Gslc(^) is concave, as shown because it represents a stable state. 

E. Constant internal order parameter 

Let us first consider the additional variable ^ to be constant for the sake of simplicity of 
the discussion. Therefore, we will suppress ^ in this section. At absolute zero, the Gibbs 
free energy in either view takes the same value Ggl(O), as noted above; see the sohd red 
curve and the dashed blue curve in Fig. HI Thus, there are two possible Gibbs free energy 
curves below the glass transition at Tg: Ggl,uv(^) or Ggl,cv(7'). Since the entropy S 
is given by Eq. (IMj) . it is evident that the magnitude of the slope of Ggl,uv(^) will be 
smaller (~ 0) than that of Ggl,cy(T) near absolute zero because of the entropy reduction, 
as shown schematically. With no entropy reduction, Ggl,cv(^) will have a large, negative 
slope equal to — at absolute zero. Hence it will drop faster than the upper red curve, 
so that Ggl.cv(^) will be given by the short dashed blue curve in Fig. HI This curve 
will eventually connect to the lower solid blue curve Gslc(^) indicated by "Supercooled 
Liquid." The combined curve (short dashed blue curve Ggl,cv(^) below Tg and solid blue 
curve GsLc(^) above Tg) remain concave at all temperatures, as shown schematically in 
the figure, as there is no inflection point. The Gibbs free energy GGL,vy(T) must also be 
continuously connected with Gslc(^) along the long dashed red curve as shown in Fig. HI 
now the connection give rise to an inflection point A. 

According to Eq. fl67j) . the Gibbs free energy must be a concave function of the tempera- 
ture so that the heat capacity remains non-negative, as the solid red curve, the short dashed 
and solid blue curves in Fig. H] are. Remember that we are talking about the total heat 
capacity of the system. It is clear that being concave, Ggl,cv(^) gives a non-negative heat 
capacity. On the other hand, this is not true of Ggl,vv{T), which contains an inflection point 
at A near Tg. It is clear that any attempt to smoothly connect the solid red piece Ggl.uvI^) 
with GsLc(^) must result in an inflection, so that it cannot remain concave everywhere. It 
must have a convex piece giving a negative heat capacity at higher temperatures, as shown 
by the arrow. Let us emphasize that to connect the two solid curves by a straight piece 
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to avoid the convex piece is contrary to physics. This will result in a discontinuity in the 
entropy and, hence, a latent heat where the connection is made. This is contrary to all 
experimental evidence. In addition, over this intermediate temperature range, the heat ca- 
pacity will turn out to be zero, which is also contrary to all experimental evidence. Thus, 
the dashed piece must not be straight for a physically realizable glass. 

F. Variable Internal order parameter 

Let us now understand the significance of Eq. fl67p when a variable internal order pa- 
rameter is present. In this case, these free energy curves will turn into surfaces defined 
over the Tq — ^ plane. In particular, the inflection point A in this flgure turns into a line 
of inflection point, as is also the case with Fig. [31 Let us consider a slice of these surfaces 
obtained by flxing ^ = ^q- As this slice will also cut the line of inflection points, the slice will 
have the middle piece in which the curvature will change sign on either side of this inflection 
point. Thus, this slice will appear similar to the curves in Fig. IHwith the middle dashed red 
curve containing an inflection point at T = Ta. At temperatures lower than Ta, the Gibbs 
free energy is concave so that Eq. fl67j) is satisfled. However, at temperatures higher than 
Ta, the Gibbs free energy becomes convex so that Eq. (1671) is not satisfled. This region then 
results in a region of negative heat capacity, and this makes the system unphysical over this 
region above Tqa. Indeed, all experimental evidence collected so far for glasses have never 
ever exhibited a negative heat capacity. 

The presence of a point of inflexion at Ta, where the heat capacity will vanish, raises 
another issue. It is obvious that Ta is a positive temperature. This gives rise to another 
unphysical aspect of UV in that the glass at this temperature can be brought in thermal 
equilibrium with any medium at any temperature. No heat exchange is possible {dQp^^ = 
Cp^^dT = 0) as the system has vanishing heat capacity. This means that at A where the heat 
capacity is zero, the temperature has no physical signiflcance for a glass; see the discussion 
above. 

In addition to these problem, the negative heat capacity is in contradiction with fl69|) . 
The negative curvature of the Gibbs free energy G'gl,uv(^) implies that the entropy must 
be a decreasing function of the temperature, in direct contradiction with the behavior of 
the entropy 5'gl,uv(^) in Fig- El The entropy derived from Ggl,uv(^) is different from the 
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entropy 5'gl,uv(^) in Fig. |3l Thus, there is an internal inconsistency, a thermodynamic 
inconsistency, within the unconventional view. The lack of thermodynamic consistency and 
the violation of Eq. (1671) rules out the unconventional view as physically relevant. Only the 
conventional view can be supported by the second law. 



VIII. CONCLUSIONS 



The formulation of the second law of thermodynamics in Eq. ([2]) presupposes the exis- 
tence and continuity of the entropy 5*0 of an isolated system under all possible conditions, 
and not only when the system is in equilibrium. This point is, however, not always appre- 
ciated. A microstate is identified by specifying some extensive mechanical quantities of the 
system such as energy, volume, etc. that we have identified here as state variables denoted 
collectively by Z(t). These variables, and therefore the microstates and their probabilities 
Pi{t) exists even if the fields like the temperature, pressure etc. may not exist for the sys- 
tem. As the entropy is a state function, the entropy Sq as identified in Eq. (H^ exists for 
an isolated body under all conditions, and we have assumed it to be a continuous function 
of Zq ant t. This assumption must be made for Eq. ([2]) to have any content. We start from 
this and show that this results in the continuity of the Gibbs free energy in Eq. (!25|) for an 
open system; see Theorem [21 The Gibbs free energy exists under all possible conditions of 
the system. In terms of the Gibbs free energy, the second law for an open body is given by 
Eq. m- 

If we now assume the open system to be in internal equilibrium, we can identify its fields 
and affinities by Eq. (fT9l) . For a body under internal equilibrium, certain results have been 
obtained in Sect. IIIIl the most important being Eq. ( 134]) and Eq. (1371) under cooling when 
we have mechanical equilibrium (P = Pq). The last observation is later used to rule out any 
theory of vitrification (via cooling) that predicts an increase of entropy during relaxation. 
The statistical concept of entropy is introduced in Sect. HVl given by Gibbs [3| for an isolated 
system under all conditions. It follows form the Fundamental Axiom in this section ^ . 
Under the assumption of quasi-independence, the entropy of an open system follows from 
this statistical formulation and is given in Eq. ( 15^ . 

The principles of non-equilibrium thermodynamics are enunciated in Sect. El which 
follow from the FUNDAMENTAL Axiom. These principles reduce to accepted principles 
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of equilibrium thermodynamics. Using these principles, we explain how the entropy of a 
system should remain continuous during component confinement at a glass transition. We 
then prove the continuity of S under two distinct situations; see (A2) and (Bl). This 
continuity of S is one of the most important results of this work and is derived under the 
most possible general condition; see Theorem El The statistical entropy merely explains how 
this happens. The next most important result is captured by Theorem |3] according to which 
only Eq. f lT3|) can be justified by the second law; Eq f lT4l) cannot be justified. The discussion 
also explains the concept of the residual entropy. The residual entropy is estimated by the 
experimentally obtained value 5'cxpt(0) of the entropy, see Eq. flMl) . by applying equilibrium 
thermodynamics to the glasse, a standard practice in the field. The experimental estimate 
'S'cxpt(O) is normally a non- negative value. We give a direct proof that the residual entropy 
must be bounded from below by this experimental value, see Eq. fl62|) . thus making the 
residual entropy at least as big as the experimental estimate. Any irreversibility will only 
raise the value of the residual entropy even higher than the experimental estimate. It is 
impossible for the residual entropy to be zero if the experimenatl estimate 5'expt(0) is non- 
zero. This conclusion justifies a non-zero S'r at G in Fig. El 

Using general arguments of stability and well-established non-equilibrium thermodynam- 
ics, we have also shown that Ggl,uv(^) is internally inconsistent in Sect. IVIIDI This 
conclusion is not based on any particular formulation of the entropy. It only uses the ther- 
modynamic concept of entropy. Our conclusion then adds another argument in support 
of the reality of the residual entropy. Any entropy reduction scenario during component 
confinement violates the second law and results in an internal inconsistency. 

We wish to thank Marty Goldstein for suggesting to consider the issue of the residual 
entropy using a statistical mechanical approach and him and S.V. Nemilov for useful corre- 
spondence. 
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